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POSITION ERROR PROPAGATION IN THE SIMPLEX STRAPDOWN NAVIGATION SYSTEM 


1.0 INTRODUCTION 


This report documents the results of an analysis of the effects of deter- 
ministic error sources on position error in the simplex strapdown naviga- 
tion system. Improving the long term accuracy of the system was addressed 
in two phases: understanding and controlling the error within the system, 

and defining methods of damping the net system error through the use of an 
external reference velocity or position. The analysis of error propagation 
is an initial step in controlling the error in an undamped system. 

Review of the flight and ground data revealed error containing the Schuler 
frequency as v;ell as non-repeatable trends. The Schuler frequency indi- 
cates the presence of error components which can be modeled by deterministic 
sources such as initial misalignment of sensor biases. Small initialization 
and instrument errors are magnified by their exciting sinusoidal error os- 
cillations with periods of 84 minutes, the Schuler period. Initial level 
alignment errors and accelerometer bias errors produce bounded sinusoidally 
varying position error components^ . The only unbounded terms are those 
Involving gyro bias and azimuth error coupled with velocity. All forms of 
Schuler-periodic position error were found to be sufficiently large, to 
require update or damping capability unless the source coefficients can be 
limited to values less than those used in this analysis for misalignment 
and gyro and accelerometer bias. 

The first-order effects of the deterministic error sources were determined 
with a simple error propagator which provided plots of error time functions 
In response to various source error values. 


nr 


Leondes, Guidance and Navigation of Aerospace Vehicles, p. 145 
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2.0 ERROR PROPAGATOR CHARACTERISTICS 


The propagator demonstrated position error propagation as a function of 

time within the simplex system operating environment constraints: the 

period of interest was on the order of two hours, and the speeds during 

the test runs were less than 100 knots. The trajectory over which the 

( 2 ) 

errors were propagated was horizontal at a constant heading and speed' ' 
specified in the input data for a period of 95 minutes. The error sources 
modeled included initial level and azimuth misalignment, accelerometer 
constant bias and scale factor error and gyroscope constant bias. 

Because of the short time period of interest, cross-coupling, earth-rate 

coupling and Coriolis-correction term errors could be ignored without 

(3l 

significantly changing the error shapes and magnitudes' . The assumption 
of no correlation between accelerometer bias and initial platform tilt was 
valid because the system was aligned either optically or off-line. That is, 
when self-alignment was employed, the alignment angles were passed to the 
navigation routines via the test operator. If self-alignment were included 
in the computer program so that navigation could commence without interrup- 
tion following alignment, the error propagation from the accelerometer 
constant bias would be cancelled by that from an initial tilt in the reference. 

The error propagation equations shown in Table 1 are derived' ' for a locally 
level system, and include the transformation matrix components where they 
are needed to relate the sensor error sources from the triad to the level 
coordinate system. 

3.0 ERROR CHARACTERISTICS 

Although the propagator was designed to yield the net error resulting from 
any set of one or more source errors, the source errors were propagated 
individually to allow visibility of the relative significance of the sources. 


^^^ibid. ; p. 143 

1 3) 

' 'Pitman, Inertial Guidance, p. 166 
(41 

' 'Ibid.; pp. 169; Kayton, Fried, Avionics Navigation Systems, p. 325 
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The test-case values for each source were chosen to be within specification 
or the bounds of experience and to be a round number so that the test results 
could be easily scaled up or down in the evaluation of a subsequent case. 

Roll and pitch initial misalignment test values were selected to be 30 arc- 
seconds, representing the observed error for an optical alignment of a rel- 
atively smooth day. The least alignment error for an optimal optical align- 
ment was approximately 8 arcseconds. 

Azimuth misalignment v/as selected to be 30 arcseconds. 

-5 

The test case value for accelerometer bias was 5 x 10 G. Bias stability 

-5 0 

was specified to be 1.5 x 10 G/ F and the temperature sensitivity of the 
bias was 1.0 x 10”^G/°F. The test value represented the effect of a 5°F 
temperature change at the instrument during the period following calibration. 
The assumption of a 5°F temperature change on the sensing element may or may 
not be representative of flight performance. Although the ambient compart- 
ment air temperature, which conditioned the sensor block, did fluctuate 
more than five degrees, sensor block temperature was not instrumented. 

Accelerometer scale factor error was tested at its specification limit of 
60 ppm. 

Gyro bias test values were set at 0.01 °/hr., the one-sigma (lo) random 
drift. A temperature change of 5*^F would have caused a bias error of .05°/ 
hr. and would have quintupled the error generated in response to the bias. 

4.0 PROPAGATED ERROR 


The error summary. Table 2, shows the, error resulting from each source at 
various true headings. The supporting plots showing the shapes of the first 
order approximations are included as Figures 1 through 32. Because the 
strapdown was mechanically constrained to vehicle azimuth, the east and 
north position errors driven by a given source varied with azimuth. The 
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root-sum square (r.s.s.), or resultant, errors at peaks and ends of runs 
were relatively independent of azimuth. Because the accelerometers and 
gyros were mounted in the triad, one hardware error source propagated 
along all three axes in the L coordinate system. Only X and Y errors were 
calculated. The r.s.s. error from a particular hardware source was inde- 
pendent of azimuth or position in the triad. Care should be exercised when 
assessing the error caused by more than one gyro bias or more than one 
accelerometer bias. The resultant position error must be determined after 
component errors are derived for the source error combination. 

The most significant error source v/as gyro constant bias. A single-instrument 
one-sigma bias error, 0.01 °/hr. produced an average r.s.s. error of approx- 
imately 4200 ft. (about 0.7 n.mi.). If the sensor temperature were five 
degrees different from its value during calibration and alignment and the 
gyro temperature sensitivity were as specified, the single-instrument bias 
contribution to resultant position error v/ould be 3.5 nautical miles. 

If all gyro bias errors were equal and if all had the same sign, position 
error would be small, qn the order of 400 feet for the cases shown. (Had the 
gyros been mounted in a way other than simplex-system triad, resultant error 
could have been quite different.) If they were equal and one had an opposite 
sign, the resultant position error contribution v;ould be increased to approx- 
imately 8500 feet (about 1.4 n.mi.), double the single-instrument contribution 
for one-sigma constant bias. This latter error would probably be increased 
by a temperature change but the assumption of quintupling is unw.irranted 
because of multiple instrument involvement. 

Accelerometer bias and initial level misalignment produce correlated 
bounded error with sinusoidal maxima at 42 minutes, the Schuler half-period. 
The average single-instrument selected accelerometer constant bias values 
was approximately 1710 feet (about 0.3 n.mi.). If temperature were not 
a factor, the specification bias stability could cause as much as 640 feet 
of resultant position error. If, as would be more representative of 
actual operation, more than one accelerometer bias error existed, the 
resultant error would be affected similarly to the way in which gyro-bias 
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resultant error varied. If all accelerometer bias were equal with the 
same sign, the resultant error would be about 3 % of the resultant produced 
by a single-instrument bias error; if they were equal but one had the sign 
opposite that of the other two, the resultant error would be double the 
single value. 

The effect of misalignment is not as complex because the triad is not 
directly involved - the propagation is somewhat more pure. The maximum 
resultant position error caused by misalignment about pitch or about roll 
was approximately 6090 feet (about 1.0 n.mi.). If roll and pitch misalign- 
ment or levelling error were both present, the effect v;ould add component- 
wise to affect the change in resultant position error. If, in that 
hypothetical case, roll error were controllable to 15 arcseconds, and pitch 
error remained about 30 arcseconds, the resultant error v/ould be about 12 % 
greater than that from a single misalignment. Opposite signs would change 
the component error but not the resultant magnitude. 

Most of the cases were initialized with the sensor-cube known off-level 
angles. Beta (b) and Theta (e), at zero. Two cases with non-zero initial 
level angles verified that the error propagation caused by gyro bias was 
affected by the initial tilt of the triad. 

Azimuth error and scale factor error are insignificant for the cases exe- 
cuted because of the relatively low speed of 80 knots. 

The plots are included as supportive and reference material. They are 
keyed to the Summary Table (Table 2) by Run Number. 

5.0 SUMMARY 


The significant deterministic error sources in the simplex strapdown navi- 
gation system are initial misalignment, gyro bias and accelerometer bias. 
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With the error coefficients selected for analysis, the system could accumu- 
late several miles of position error in one hour. The coefficients during 
the flight tests could have been larger than those selected. 

Self-alignment with instantaneous transfer to the navigation mode can 
reliably reduce initial misalignments to the range of 4-7 arcseconds with a 
subsequent reduction by an order of magnitude in the error propagated. In 
addition, the level errors v/ould then tend to be correlated with acceler- 
ometer bias, so that one cancelled the other. 

Gyro bias error is the most significant contributor of those identified. 

The effect of temperature change on gyro bias is known. The actual tem- 
perature value and variations for the test runs arc not known. It has 
been shown that a modest change of the sensor temperature can produce 
significant error. Temperature should, at least, be monitored during de- 
velopmental system test. Laboratory tests should also be devised to pre- 
cisely map the reaction of gyros to temperature. If significant reaction 
is determined, and significant temperature changes are detected in system 
test, software compensation for real-time temperature change should be 
provided. 

The position error propagated from gyro bias equivalent to the one-sigma 
random drift value cannot be compensated within the undamped system. Addi- 
tionally, the presence of stationary zero-mean noise in the gyro signal 

(51 

causes unbounded system error' . The system must be damped with an external 
reference velocity or position to provide satisfactory long-term accuracy. 

The results of this deterministic-error propagation can be used in refining 
the undamped system to simplify the damping requirement and process. 




Leondes, p. 146 
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Table 1 . Position Error Propagation Equations 

Initial Level Misalignment 

aX = Rq (ai|/y sin Az + cos Az) (1 - cos u^t) 

aY = Rq CCS Az - sin Az) (1 - cos u^t) 

Initial Azimuth Misalignment (ao^^ 
aX = A. ^ Vt cos Az 
aY = Ai^2 ^ 

Accelerometer Bias Error (aB^^, ABy, aB^) 

(1 - cos W 5 t) 

aX = (i.B^ + ABy ^"'SL31^ 2 

0 ,$ 

(1 - cos Ujt) 

aY = (aBj^ ^’'sL22 ^ ^SL32^ 2 

Uc 


Accelerometer Scale Factor Error (aKj^» AKy, aK^) 
aX = - ^aKj^ ^^SLll^ ^ ^SLll ^SL 12 ^ 

+ AKy (M 5 L 2 i^ V sin Az + M^^ 2 i ^SL 22 ^ 

+ aK^ (Msl 3 i^ V sin Az + M^j _32 ^ Azj 



sin ui^t 
“ 
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Table 1. Position Error Propagation Equations (Continued) 


r 


c 


4Y . - [iKj; V s1n Az + M5,_j22 V cos Az) 

+ AKy ^sL 22 sin Az + \n2 ^ 

+ AK^ (M5L31 f'^sL32 ^ ■*■ '"'sL32 


t - 


sin ujt 


Gyro Bias Error A^iyi A412) 


Level Components 

I sin (u<.t\ 

AX = Rq ^^SL22 * ^SL32^ ^ f 

/ sin u<.t\ 

"SLII * ”SL21 * - -^ 7 ^) 


Azimuth Error Components 


AX = V sin Az (AO^ ^SL13 '*’ ^^3L23 ^SL33^ 


2 2 (1 - cos u^t) 


« 2 (1 - cos ^f.t) 

AY = V cos Az (A«^ Mjli 3 * AO^ * AC^ ( t“= — 




Term Definitions 


( 


“S 

V 

Az 



earth radius 

Schuler frequency (rad.) 

Ground speed 
Azimuth 

MSL matrix (as defined in Program Requirements Ooc.) 
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